APIT/OAE/AA/81D-30 


r 


AN,  EXPERIMENTAL  INVESTIGATION  OP  THE 
PLOW  FIELD  OF  AN  EJECTOR  WING  DESIGN 
EMPLOYING  A  PHOTON  CORRELATION 
LASER  VELOCIMETER 


THESIS 


APIT/OAE/AA/81D-30 


DAMON  G.  STEPHENS 
Captain  USAF 


DT1C 

SELECTS' 

FED  1  9  1982 

A 


Approved  for  public  ralaaaoi  distribution  unlimited 


[Pll  Redacted] 


AFIT/G AE/AA/8 1 D- 30 


AN  EXPERIMENTAL  INVESTIGATION  OF  THE 
FLOW  FIELD  OF  AN  EJECTOR  WING  DESIGN 
EMPLOYING  A  PHOTON  CORRELATION  LASER  VELOCIMETER 

THESIS 

Presented  to  the  Faculty  of  the  School  of  Engineering 
of  the  Air  Force  Institute  of  Technology 
Air  University 

in  Partial  Fulfillment  of  the  Requirements 
for  the  Degree  of  Master  of  Science 

by 

Damon  G.  Stephens,  B.S. 

Captain  USAF 

Graduate  Aeronautical  Engineering 
December  1981 

\  <***'-'' 

Approved  for  public  release;  distribution  unlimited  - 


PREFACE 


This  tnesis  arose  out  of  a  need  by  the  Air  Force  Flignt 
Dynamics  Laboratory  (AFFDL)  for  mean  velocity  and  turbulence 
intensity  data  in  the  flow  field  of  an  ejector  wing  design. 

The  design  was  conceived  by  the  Vought  Corporation  while 
under  contract  with  AFFDL  and  represents  a  unique  and  exciting 
airfoil/ejector  mating 

As  a  novice  to  the  world  of  researcn  and  technical  writ¬ 
ing,  I  greatly  appreciated  the  guidance  and  expert  assistance 
of  many  talented  individuals.  My  sincere  thanks  to  Dr.  Harold 
Wright,  my  principal  faculty  advisor,  for  the  opportunity  to 
work  on  this  project  and  for  his  assistance  in  dealing  with  the 
many  equipment  problems  that  impeded  its  progress.  Sincere 
thanks  also  to  Dr.  George  Catalano  for  sharing  with  me  his 
wealth  of  laser  velocimeter  experimental  experience.  I  also 
wish  to  thank  Maj.  Mike  Smith  for  his  assistance  in  finding 
answers  to  questions  that  arose. 

For  their  efforts  in  setting  up  and  maintaining  the  ex¬ 
perimental  apparatus  I  wish  to  thank  Mr.  William  Eaker  and 
Mr.  Leroy  Cannon.  Also  a  word  of  appreciation  and  thanks  to 
Mr.  Nick  Yardich  and  Mr.  Scotty  Whitt  for  their  dedication  to 
keeping  the  "Blue  Beast"  (AFIT  Smoke  Tunnel)  in  operation 
throughout  my  research.  And  for  their  useful  suggestions  and 
fine  craf traanship  in  the  construction  and  modification  of  the 
ejector  wing  model,  I  wish  to  thank  Mr.  Carl  Shortt,  Mr.  Russel 
Murray,  and  Mr.  Jack  Tiffany  of  the  AFIT  School  Shop.  Also, 
for  their  excellent  photographic  support  I  wish  to  thank  Mr. 


ii 


Ed  Fields  and  Mr.  Dave  Cunningham  of  the  Technical  Photo 
Division. 

In  addition,  and  most  importantly,  I  wish  to  express  my 
deepest  gratitude  to  my  wife,  Kathy,  for  her  patience,  love, 
and  understanding  during  the  course  of  study  required  in  the 
preparation  of  this  thesis. 


iii 


CONTENTS 


”  Page 

Preface . . . . .  ii 

List  of  Figures . .  v 

List  of  Symbols . vii 

Abstract . . .  ix 

I.  Introduction . . .  1 

Background .  1 

High  Lift  Configurations . .  2 

Objectives.... .  3 

Approach. .  3 

II.  Experimental  Apparatus  and  Procedure .  5 

Flow  System............ .  5 

Ejector  Wing  Design  Model .  7 

Laser  Doppler  Velocimeter .  14 

III.  Sources  or  Error  and  Equipment  Limitations .  20 

IV.  Results  and  Analysis . . .  23 

Flow  Visualization .  23 

Laser  Velocimeter...,. .  32 

AO  A  -  15° .  32 

AOA  =  0° .  64 

V.  Conclusions . . . . .  60 

VI.  Recommendations .  63 

Bibliography..., . . .  86 

Appendix  A:  Autocorrelation  Function .  88 

Appendix  B:  Mean  Velocity .  94 

Appendix  C:  Turbulence  Intensity . .  98 

Appendix  D:  Sample  Calculations .  101 

Vita .  109 


LIST  OP  FIGURES 


Figure  Title  Page 

1  Flow  System  6 

2  Ejector  Wing  Design  b 

3  Mercury  Manometer  Calibration 

Setup  9 

k  Ejector  11 

5  Ejector  Velocity  Profile  12 

6  Ejector  Velocity/Mercury  Manometer 

Calibration  12 

7  Laser  Velocimeter  System  (Top  View)  15 

b  Optical  Bench  Showing  Optical 

Path  (Side  View)  15 

9  Laser  Beam  Control  Volume  16 

10  Oscilloscope  Autocorrelation 

Display  17 

11  Pictorial  View  of  Experimental 

Apparatus  19 

12  Laser  Beams  Oblique  to  Airfoil  20 

13-26  Flow  Visualization  Photographs  25-31 

27  Leading  Edge  Stagnation  Point 

Movement  3b 

2b  Survey  Locations,  AOA  =15°  37 

29-1*1  Mean  Velocity  Profiles,  AOA  =  15°  3&-50 

42- 54  Turbulence  Intensity  Profiles, 

AOA  =  15°  51-63 

55  Wake  Survey  66 

56  Survey  Locations,  AOA  =  0°  6b 

57-62  Mean  Velocity  Profiles,  AOA  =  0°  69-7L 


v 


Title 


Turbulence  Intensity  Profiles, 

AO A  =  0°  75-79 

Autocorrelation  of  u(t)  and  u(t+*t  )  8b 

Autocorrelation  Function  89 

Oscilloscope  Displayed  Autocorrela¬ 
tion  Function  91 

Displayed  Skewed  Autocorrelation 
Function  92 

Replotted  Unskewed  Autocorrelation 
Function  92 

Autocorrelation  Function  93 

Mean  Velocity  914 

Illustration  of  L,  D,  e  95 

Autocorrelation  Peaks  and  Valleys  96 

Velocity  Fluctuations  9b 

High  Turbulence  Displayed 
Autocorrelation  Function  102 

Plot  of  Displayed  Skewed 
Autocorrelation  Function (Sample 
Calculations)  103 

Least  Squares  Fit(Sample  Calculations)  101+ 

Least  Squares  Fit(Sample  Calculations)  106 

Unskewed  Autocorrelation  Function 
(Sample  Calculations)  107 

Nondimensional  Autocorrelation 
Function(Sample  Calculations)  107 

vi 


SYMBOLS 


« 


Symbol 


Description 


Units 


AOA 

C 

D 


D' 


ETA 

*d 


angle  of  attack 

degrees 

chord  length  of  model 

cm 

distance  the  two  laser  beams 
are  apart  at  a  perpendicular 
reference  surface (See  L) 

inches 

drag  per  unit  span 

lbf/ft 

turbulence  Intensity 

doppler  shift  frequency 

KHz 

doppler  shift  frequency 
with  Phase  Modulation 

Khz 

Af 

frequency  shift 

imposed  on 

the  two  laser  beams  by  the 

Phase  Modulator 

G(t) 

autocorrelation 

function 

G1 

first  valley  in 
function 

the  autocorrelation 

G2  first  peak  after  the  first  valley 

(G1 )  in  the  autocorrelation  function 

G3  second  valley  in  the  autocorrelation 

function 


L 


LDV 

N 


distance  from  the  laser  beam  inches 

intersection  point  (control 

volume)  to  a  perpendicular 

reference  surface  upon  which 

the  beams  shine 

laser  doppler  velocimeter 

number  of  fringes  in  half  of 
control  volume 


PM  Tube 

♦ 


photomultiplier  tube 


vii 


Symbol 


Description 


ro 

S 

T 

Ue 


laser  beam  radius 
fringe  spacing 
sample  time 
ejector  velocity 


tVUfs 


ejector  blowing  ratio 


f 

mX 


URMS 

X 


x' 

y' 

z ' 


free  stream  velocity 

mean  velocity  component 
parallel  to  the  free 
stream 

mean  veloctiy  component 
parallel  to  the  free 
stream  with  Phase  Modulation 

root  mean  square  velocity 

direction  parallel  to  free 
stream 

chordwise  length  dimension 
from  leading  edge  of  airfoil 

spanwise  length  dimension 

vertical  length  dimension 
from  the  chordline  of  airfoil 


z  vertical  length  dimension 

measured  from  the  surface  of 
airfoil  at  each  survey  location 

9  half  angle  between  the  inter¬ 

secting  laser  beams 

X  wavelength  of  the  he-ne  laser 


U  refractive  Index  of  Air 

f  air  density 


T 


delay  time 


Units 

mm 

m 

sec 

m/sec 

m/sec 

m/sec 

m/sec 

m/sec 

cm 

cm 

cm 

cm 

degrees 

m 

lbm/ft3 

sec 


viii 


AFIT/GAE/ AA/8l D-30 


ABSTRACT 

The  flow  field  about  a  dual  element  airfoil  model  em¬ 
ploying  an  ejector  for  aerodynamic  blowing  was  investigated. 
Flow  visualization  was  obtained  by  smoke  tunnel  testing. 

Mean  velocities  (mean  velocity  components  parallel  to  the 
free  stream),  and  turbulence  intensities  were  determined 
at  various  flow  field  locations.  All  data  was  obtained 
through  the  use  of  a  laser  doppler  velocimeter  (LDV)  using 
a  photon  correlation  processing  scheme.  Flow  field  proper¬ 
ties  were  computed  from  the  LDV  generated  autocorrelation 
function. 

The  free  stream  velocity  and  Reynolds  number,  based  on 
model  chord  lengtn,  were  8  m/sec.  and  325,000,  respectively. 

The  ejector  wing  model  was  tested  with  Ue jector/Ufree  stream 
(U  /U  )  =  2 .0  and  no  ejector  flow.  The  U  /TJ  =2.0  blowing 

0X3  0X3 

ratio  delivers  0.046  lb  /sec  or  0.60  ft'Vsec  at  the  ejector 

m 

face.  Flow  visualization  photography  was  conducted  at  seven 
angles  of  attack:  -5°,  0°,  5°,  10°,  15°,  20°,  and  25°  and 
LDV  data  was  acquired  at  0°  and  15°  angles  of  attack. 

The  results  are  presented  in  the  form  of  flow  visual¬ 
ization  photographs  and  profiles  of  mean  velocities  and  turb¬ 
ulence  intensities.  The  ejector  wing  design  achieves  superior 
performance  (enhanced  lift,  reduced  drag)  over  an  equivalent 
solid  airfoil. 


I.  INTRODUCTION 


Background 

The  ejector  wing  design  analyzed  in  this  paper  was 
conceived  by  the  Vought  Corporation  while  under  contract 
with  the  Air  Force  Flight  Dynamics  Laboratory  (AFFDL)  at 
Wright-Patteraon  AFB,  Ohio  (Reference  11).  The  model  was 
fabricated  and  mounted  in  the  Air  Force  Institute  of  Tech¬ 
nology  (AFIT)  two-dimensional  smoke  tunnel  at  Wright-Pat- 
terson  AFB  prior  to  this  present  investigation.  A  sche¬ 
matic  of  the  ejector  wing  design  is  presented  in  Figure  2. 
The  thrust  of  this  project  is  to  study  the  flow  field  about 
the  model  as  an  aid  to  the  Flight  Dynamics  Laboratory's  in¬ 
vestigation  of  the  ejector  wing  design  for  possible  V/STCL 
application. 

The  Air  Force  Institute  of  Technology's  Laser  Doppler 
Velocimeter  (LDV)  system  incorporating  a  photon  correlation 
processing  scheme  is  used  to  obtain  the  experimental  data. 
Flow  visualization  through  kerosene  smoke  generation  is 
used  to  assist  in  determining  the  overall  flow  field  var¬ 
iations  caused  by  ejector  operation.  There  is  a  great  deal 
of  theoretical  analysis  and  experimental  documentation  of 
simple  ejectors  available  in  the  literature;  however,  very 
little  information  is  available  concerning  the  flow  field 
about  a  wing  incorporating  an  ejector.  Hopefully,  an  ex¬ 
perimental  investigation  of  this  unique  design  will  contrib¬ 
ute  to  a  better  understanding  of  this  concept. 
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High  Lift  Configurations 


Multi-element  airfoils  have  long  been  recognized  as 
beneficial  in  obtaining  increased  lift  and  reduced  drag. 
The  benefits  include:  increased  circulation,  decreased 
wake  region,  and  enhanced  boundary  layer  adhesion.  The 
increased  circulation  is  due  to  the  location  of  the  down¬ 
stream  element  near  tne  rear  of  the  upstream  element.  The 
constriction  to  flow  and  the  circulation  about  tne  down¬ 
stream  element  places  the  trailing  edge  of  the  upstream 
element  in  a  region  of  high  velocity.  In  order  for  the 
flow  on  the  upper  surface  of  the  upstream  element  to  meet 
the  decreased  pressure  at  its  trailing  edge  it  must  ac¬ 
celerate.  Tnis  serves  to  increase  the  circulation  about 
the  entire  airfoil  and  decrease  the  adverse  pressure  grad¬ 
ient  on  the  upper  surface  of  the  upstream  element.  Thus, 
lift  is  increased  and,  since  separation  is  delayed,  the 
wake  region  and  its  associated  drag  are  decreased.  The 
introduction  of  an  ejector  blowing  high  velocity  air  into 
the  constricted  duct  and  over  the  upper  surface  of  the 
downstream  element  further  enhances  these  effects.  Also, 
the  high  energy  air  from  the  ejector  aids  in  keeping  the 
downstream  element’s  boundary  layer  attached.  For  further 
information  on  high  lift  devices,  an  excellent  article  by 
A.O.M.  Smith  (Reference  19)  should  be  consulted. 
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The  objectives  of  this  study  are: 


1 .  Using  the  two  dimensional  AFIT  smoke  tunnel  obtain  flow 
visualization  photographs  of  the  ejector  wing  model  at  var¬ 
ious  angles  of  attack:  -5°,  0°,  5°,  10°,  15°,  20°,  25°,  and 
ejector  blowing  ratios:  Ue/U^.g  =  0.0  and  2.0. 

2.  Using  a  Laser  Doppler  Velocimeter  obtain  and  plot  mean 
velocity  and  turbulence  intensity  profiles  at  various  flow 
field  locations. 

3.  From  the  acquired  data  draw  conclusions  as  to  the  aero¬ 
dynamic  benefit  of  this  ejector  wing  design. 

Approach 

In  order  to  meet  the  objectives  of  this  study  two  dif¬ 
ferent  investigative  schemes  will  be  conducted: 

1 .  Flow  visualization  photography  through  the  ubc  of 
kerosene  smoke. 

2.  Mean  velocity  and  turbulence  intensity  data  at  de¬ 
signated  flow  field  locations  obtained  by  a  LDV. 

Both  investigative  schemes  will  be  employed  with  smoke 

tunnel  free  stream  velocity  (Ufs)  of  6  m/sec  and  ejector 

blowing  ratios,  Ue/Ufa»  of  0.0  and  2.0.  Flow  visualization 

photographs  will  be  taken  at  -5°,  0°,  5°,  10°,  15°,  20°,  and 

25°  angles  of  attack. 

LDV  data  will  be  obtained  at  0°  and  15°  angles  of  attack. 
At  15°  angle  of  attack  data  will  be  taken  with  the  ejector 
duct  open.  At  0°  angle  of  attack,  a  comparison  of  the  ejector 


wing  design  with  an  equivalent  solid  airfoil  will  be  con¬ 


ducted.  Also,  mean  velocity  wake  surveys  will  be  obtained 
and  plotted  in  an  attempt  to  determine  relative  drag. 


k 


II.  EXPERIMENTAL  APPARATUS  AND  PROCEDURE 


Flow  System 

In  order  to  provide  a  uniform  free  stream  flow,  the  Air 
Force  Institute  of  Technology's  two  dimensional  smoke  tunnel 
(Figure  1 )  is  used.  The  smoke  tunnel  also  provides  the  flow 
field  visualization  capability  as  well  as  an  artifically  seeded 
flow  compatable  with  the  LDV  system.  Maximum  tunnel  velocity 
is  23  m/sec  using  an  open  return  flow  system  obtained  by  two 
diffuser  isolated  1.5  hp  motors.  For  this  study  a  free  stream 
velocity  of  8  m/sec  and  a  Reynolds  Number  based  on  model  chord 
length  of  325,000  is  established.  This  velocity  is  obtained 
and  monitored  by  a  Prandtl  type  pitot-static  probe  and  a  20  inch 
micromanometer.  The  free  stream  velocity  is  maintainable  to 
within  +  0.1  m/sec  or  within  1.25%  of  the  established  free 
stream  velocity.  The  test  section  measures  1 .5  m  in  length, 
1.0m  in  height,  and  .07  m  in  width.  The  front  wall  is  remov¬ 
able  for  ease  of  model  access  and  its  frame  can  accomodate  any 
0.635  cm  thick  glass  material,  such  as,  plexiglass  or  plate- 
glass.  The  rear  wall,  into  which  the  model  is  mounted,  is 
0.635  cm  laminated  plateglass.  Smoke  (kerosene  vapor)  is  gen¬ 
erated  by  two  900  watt  inconel  heaters  which  boil  the  kerosene 
fuel  at  910°R.  The  vapor  particles  are  then  mixed  with  cool 
compressed  laboratory  air,  passed  through  a  water  condenser, 
and  injected  into  the  flow  field  upstream  of  the  test  section. 
The  dense,  white  smoke  produced  is  nontoxic,  noncorrosive,  and 
provides  excellent  visualization  capabilities. 
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Condensation 

Chamber 


Ejector  Wing  Design  Model 

The  model  consists  of  a  dual  element  airfoil  with  an 
incorporated  ejector  (Figure  2  ).  It  is  constructed  of  wood 
and  aluminum  with  a  smooth  flat  black  surface  to  minimize 
light  scattering.  The  forward  and  aft  airfoils  and  ejector 
are  permanently  secured  to  a  plexiglass  mounting  plate  and 
thus  remain  in  the  same  configuration  relative  to  each  other 
as  shown.  The  free  space  area  or  duct  between  the  forward 
and  aft  airfoils  where  the  ejector  is  located  is  referred  to 
as  the  ejector  duct.  The  entire  airfoil  is  mounted  to  tne 
back  wall  of  the  smoke  tunnel  at  approximately  the  center  of 
the  test  section.  Two  dimensional  flow  is  assured  by  closing 
the  test  section  securely  against  the  model  about  its  mounting 
apparatus.  Test  section  blockage  caused  by  the  model  ranges 
from  11%  at  0°  angle  of  attack  to  29%  at  25°  angle  of  attack. 
The  air  is  blown  into  the  ejector  in  a  spanwise  direction  and 
ejected  in  a  chordwise  direction  at  the  ejector  face.  The 
laboratory  compressed  air  reservior  supplies  the  ejector  en¬ 
abling  various  flow  rates  to  be  achieved. 

The  first  task  in  this  study  was  to  determine  the  ejector 
velocity  (Ue).  A  pitot-static  tube  was  inserted  in  the  front 
wall  of  the  test  section  to  probe  the  center  of  the  ejector 
face  in  a  spanwise  direction  (Figure  3  ).  This  probe  was  con¬ 
ducted  with  the  model  at  0°  angle  of  attack.  Since  the  pitot 


tube  would  interfere  with  the  flow  field  and  laser  data  acqui¬ 
sition,  a  pressure  tap  connected  to  a  mercury  manometer  was 


Figure  2  Ejector  Wing  Design 


Figure  3  Mercury  Manometer  Calibration  >etup 
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placed  in  the  blowing  line  (air  line).  By  determining 

the  ejector  velocity  at  various  pressures,  the  manometer 

could  be  calibrated  and  its  readings  used  to  reestablish 

desired  ejector  velocities  during  the  investigation.  In 

addition  to  determining  the  ejector  velocity,  the  pitot 

tube  setup  was  used  to  modify  the  original  ejector  design 

to  achieve  a  nearly  uniform  spanwise  velocity  profile  at 

the  ejector  face.  The  new  design  was  needed  due  to  the 

nonuniform  characteristics  of  right  angle  blowing,  i.e., 

excessive  ejector  face  velocities  at  the  outboard  section 

and  near  stagnation  conditions  at  the  inboard  section.  Use 

of  a  blowing  tube  with  decreasing  diameter  holes  inboard 

to  outboard  with  a  screen  downstream  of  the  tube  provided 

the  desired  uniform  spanwise  velocity  profile  (Figure  I*.), 

Blockage  of  the  ejector  face  caused  by  the  screen  is  67%, 

allowing  back  pressure  to  develop  which  aided  in  yielding 

the  nearly  uniform  velocity  profile.  Figure  5  shows  the 

velocity  profile  for  U  =  16.0  m/sec  which  is  the  ejector 

© 

on  velocity  used  throughout  this  investigation  = 

2.0).  After  achieving  the  desired  velocity  profile,  the 
mercury  manometer  was  then  calibrated  (Figure  6). 

In  addition  to  ejector  internal  modifications,  a 
method  for  changing  the  angle  of  attack  of  the  model  from 
outside  the  test  section  was  devised  and  implemented. 
Rotation  of  the  model  exterior  to  the  test  section  was 
accomplished  by  the  use  of  strong  magnets.  Several  strong 
magnets  were  embedded  in  the  aft  airfoil  flush  with  the 
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Figure  U  Ejector 


plexiglass  mounting  plate.  Additional  magnets  were  then 
used  to  rotate  the  model  from  outside  of  the  test  section. 
By  drawing  reference  lines  on  the  back  wall  of  the  test 
section,  various  angles  of  attack  could  be  achieved. 
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Laser  Doppler  Veloclmeter 

In  order  to  obtain  the  mean  velocity  and  turbulence 
intensity  profiles  for  this  investigation  a  Malvern  Doppler 
Velocimeter  (LDV)  was  employed  (Figures  ?  and  8).  The  ex¬ 
perimental  benefit  of  the  LDV  system  for  determination  of 
flow  properties  is  due  to  its  inherent  noninterference  with 
the  flow.  In  addition,  no  calibration  is  required.  The 
laser  beams,  unlike  hot  wire  anemometers,  pitot  tubes,  and 
pressure  transducers  do  not  create  disturbances  or  pressure 
losses  and  do  not  suffer  damage  from  the  flow  enviornment. 
However,  the  correct  optical  alignment  is  critical  and  thus 
requires  constant  monitoring.  Also,  light  paths  into  and 
out  of  a  transparent  fluid  containing  suitable  light  scat¬ 
tering  particles  is  required. 

The  Laser  Velocimeter  consisting  of  a  Spectra-Physics 
15  milliwatt  helium-neon  laser,  Malvern  beamsplitter  and 
phase  modulator,  two  silvered  plane  mirrors,  a  100  cm  focal 
length  convex  focusing  lens  ,  Malvern  photomultiplier  tube 
(PM  Tube),  Malvern  digital  photon  correlator,  oscilloscope, 
supporting  optical  benches,  and  associated  power  units  was 
constructed  for  a  previous  AFIT  master's  thesis  by  Walterick 
(Reference  23). 

A  brief  outline  of  the  velocimeter' s  operation  is  pre¬ 
sented.  For  more  detailed  operation  and  theory  References 
16,  17,  and  23  should  be  consulted. 

The  he-ne  laser  provides  a  single  1.1  mm  diameter  mono¬ 
chromatic  (  A  =  6328yttm)  red  beam  which  is  split  into  two 
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Figure  7  Laser  Velocimeter  System  (Top  View) 


Figure  8  Optical  Bench  Showing  Optical  Path  (Side  View) 
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beams  by  a  beamsplitter.  The  two  beams  pass  through  a 
phase  modulator  and  are  then  transmitted  along  the  optical 
path  to  the  convex  focusing  lens  by  two  silvered  plane 
mirrors.  As  a  result  of  passing  through  the  lens  the 
beams  intersect  at  the  focal  point  and  form  an  ellipsoidal 
control  volume  composed  of  alternating  constructive  (light) 
and  destructive  (dark)  interference  fringes  (Figure  9  ). 


Mean  velocity  and  turbulence  intensity  data  can  then  be 
obtained  by  focusing  a  photomultiplier  tube  (PM  tube)  on  the 
control  volume  as  naturally  occuring  particles  augmented  by 
kerosene  vapor  pass  through  these  light  and  dark  fringes. 

Since  the  control  volume  fringes  are  oriented  perpendicular 
to  the  free  stream  flow,  only  mean  velocity  components  parallel 
to  the  free  stream  can  be  obtained.  The  PM  tube  detects  the 


16 


the  photons  scattered  by  the  particles  and  in  conjunction 
with  a  digital  photon  correlator  generates  an  autocorrelation 
function  which  is  displayed  on  an  oscilloscope  (Figure  10). 
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Figure  10  Oscilloscope  Autocorrelation 
Display 

In  order  to  generate  an  adequate  autocorrelation  function, 
i.e.,  one  with  a  damped  sinusoidal  shape,  under  certain  adverse 
conditions  a  phase  modulator  unit  is  employed.  Flow  which  is 
too  slow,  too  fast,  or  too  turbulent  cannot  be  adequately  ev¬ 
aluated  without  the  phase  modulator.  Also,  flow  direction 
cannot  be  determined  unless  a  phase  modulator  is  used.  V'ith- 
out  phase  modulation  the  interference  fringes  in  the  control 
volume  are  stationary  with  respect  to  the  flow.  By  employing 
a  phase  modulator  the  fringes  can  be  made  to  move  either  with 
or  against  the  flow  by  applying  an  appropriate  frequency  shift 
to  each  beam.  Attainable  frequency  shifts  are  20  KHz  to  1  MHz. 
By  moving  the  fringes  with  the  flow,  for  example,  a  fast  flow 
can  be  "slowed"  allowing  more  photon  scattering  and  thus 
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improving  the  PM  tube's  efficiency  and  yielding  a  better 
autocorrelation  function.  In  addition,  by  controlling 
the  direction  of  fringe  movement  flow  direction  can  be 
evaluated. 

All  data  in  this  investigation  was  taken  at  the 
spanwise  (y-direction)  midplane  of  the  test  section  (see 
Figure  7).  Therefore,  the  LDV  control  volume  was  positioned 
at  the  test  section  midplane.  Traversing  of  the  control 
volume  in  the  chordwise  (x-direction)  is  accomplished  by 
translation  of  the  entire  bench.  Traverses  in  the  z- 
direction  are  accomplished  by  translation  of  a  secondary 
table  through  a  chain  and  sprocket  mechanism  (Figures  7 
and  6 ) . 

A  brief  explanation  and  discussion  of  the  autocorrel¬ 
ation  function  plus  definitions  and  discussions  of  turbulence 
properties  are  presented  in  the  Appendix.  In  addition,  the 
reduction  of  autocorrelation  data,  both  with  and  without 
phase  modulator  employment,  into  mean  velocities  and  turb¬ 
ulence  intensities  is  presented  in  Appendices  B  and  C. 

A  pictorial  view  of  the  entire  experimental  apparatus 
is  presented  in  Figure  1 1 . 
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III.  SOURCES  OP  ERROR  AND  EQUIPMENT  LIMITATIONS 

Before  presenting  the  results  and  subsequent  conclu¬ 
sions  and  recommendations,  the  possible  sources  of  error 
and  limitations  of  the  equipment  are  presented  below: 

1 .  The  LDV  optical  alignment,  in  the  absence  of  a 
rigid  optical  table,  was  very  sensitive  and  required 
constant  monitoring  and  frequent  realignment. 

2.  Background  light  which  results  in  optical  noise, 
although  minimized  to  the  greatest  extent  possible, 
may  still  have  yielded  excessive  noise  and  adversely 
affected  the  data. 

3.  It  was  found  that  data  could  not  be  obtained  on 
the  surface  of  the  model  nor  within  2-3  mm  of  the 
surface  due  to  the  following: 

a.  The  path  at  which  the  laser  beams  enter  the 
teat  section  is  not  parallel  to  the  airfoil's 
surface  as  a  result  of  floor/optical  bench  ir¬ 
regularities  (Figure  12). 


Figure  12  Laser  Beams  Oblique  to  the  Airfoil 
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If  the  beams  are  lowered  below  that  shown  in 
Figure  12,  the  airfoil  surface  would  block  them. 

b.  Near  the  surface  the  beams  caused  excessive 
reflections  which  tended  to  overpower  the  PM 
tube  and  mask  usable  data. 

c.  Occasionally,  the  turbulence  at  the  surface 
was  in  excess  of  the  system's  capability.  This 
indicated  a  need  for  higher  phase  modulation  than 
1  MHz. 

As  a  result  of  these  limitations,  data  nearer  than  2 
to  3  mi  from  the  surface  was  impossible  to  obtain. 
Therefore,  useful  data  within  the  boundary  layer  could 
not  be  accurately  acquired.  In  addition,  all  mean  vel¬ 
ocities  obtained  are  not  total  velocity  vectors  but 
only  components  aligned  with  the  free  stream.  There¬ 
fore,  no  meaningful  boundary  layer  analysis  could  be 
conducted  using  this  equipment  and  this  model  size. 

4.  Due  to  limitations  imposed  by  the  optical  benches, 
the  maximum  downward  (negative  z-direction)  movement 
of  the  control  volume  was  approximately  12  cm  below 
the  center  of  the  test  section. 

5.  Initially,  plexiglass,  used  for  prior  research  ap¬ 
plications  for  safety  reasons,  was  used  as  the  front 
wall  of  the  smoke  tunnel's  test  section.  However,  due 
to  the  softness  of  the  plexiglass  material  it  was  easily 
scratched  during  normal  use.  These  scratch  marks  along 
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with  the  refractive  variations  of  the  material  caused 
by  stresses  induced  a  large  scattering  of  laser  light 
which,  at  times,  generated  overwhelming  optical  noise. 
Therefore,  the  plexiglass  was  replaced  by  plate  glass. 
This  reduced  the  optical  noise  and  yielded  a  better 
signal. 
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IV.  RESULTS  AND  ANALYSIS 


The  results  are  presented  in  the  order  in  which  they 
were  experimentally  obtained: 

1 .  Plow  Visualization  Photographs 

2.  LDV  Experimental  Data 

The  main  purpose  of  obtaining  flow  visualization  photo¬ 
graphs  is  to  gain  an  appreciation  of  the  overall  changes  in 
the  flow  field  due  to  ejector  employment.  In  addition,  the 
photographs  serve  as  a  very  useful  aid  in  determining  the 
locations  at  which  LDV  data  should  be  acquired.  Hence,  the 
reason  for  the  experimental  sequence. 

Flow  Visualization 

The  visual  effects  of  ejector  operation  on  the  flow 
field  are  documented  in  Figures  13  through  26.  In  all 
photographs  the  ejector  duct  is  open.  The  arrows  in  each 
figure  (photograph)  identify  a  common  streak  line  which 
can  be  used  as  an  aid  in  flow  field  comparisons  between 
the  ejector  off  and  ejector  on  operational  modes. 

Two  observations  can  be  readily  made.  First,  the 
fluid,  here  visualized  by  streak  lines,  is  drawn  up  from 
its  original  position  before  ejector  employment.  Second, 
the  streak  lines  in  the  airfoil  wake  region  directly  above 
the  aft  airfoil  appear  less  diffused  and  more  well  defined 
with  ejector  operation.  An  excellent  illustration  of  the 
first  observation  occurs  at  AOA  =10°  (Figures  19  and  20). 
The  arrow  in  Figure  19  indicates  a  streak  line  which  flows 
under  the  forward  airfoil  and  into  the  ejector  duct  with¬ 
out  ejector  operation.  The  same  streak  line,  identified 
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by  the  arrow  in  Figure  20,  flows  over  the  upper  surface  of 
the  forward  airfoil  with  ejector  operation.  This  phenomenon, 
as  discussed  earlier,  is  attributable  to  the  ejector's  crea¬ 
tion  of  a  high  velocity  region  at  the  trailing  edge  of  the 
forward  airfoil.  This  causes  an  acceleration  of  the  flow  on 
the  upper  surface  and  a  drop  in  static  pressure  or  increase 
in  suction.  This  effect  draws  fluid  up  from  its  original 
position  before  ejector  employment.  The  net  result  is  a 
greater  pressure  differential  and  an  increase  in  effective 
angle  of  attack  yielding  greater  lift.  The  second  obser¬ 
vation  would  seem  to  indicate  a  smaller,  more  well  behaved 
wake  region  which  could  lead  to  a  drag  reduction.  An  ex¬ 
cellent  illustration  of  this  wake  region  phenomenon  occurs 
at  AOA  =  15°  (Figures  21  and  22).  Note  the  apparent  decrease 
in  wake  region  size  above  the  aft  airfoil  as  well  as  the  less 
diffused,  more  definable  flow  pattern. 

It  can  be  seen  that  the  leading  edge  stagnation  point  ap¬ 
pears  to  move  down  or  in  a  counterclockwise  direction  at  the 
nose  with  the  ejector  on.  This  iridicatea  an  increase  in  cir¬ 
culation.  A  final  observation  is  the  apperance  that  more  flow 
is  being  entrained  into  the  ejector  duct  region  during  ejector 
operation.  Note  this  effect  at  AOA  *  5°(  Figures  17  and  1 8 ) . 
More  fluid  can  be  seen  within  the  ejector  duct  with  the  ejector 
on  than  with  it  off.  This  is  an  expected  ejector  benefit  and 
should  contribute  to  enhanced  performance. 


I*  AOA  =  0°  ,  Ejector 


OFF 


Figure 
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Figure  1 7  AOA  =  5°  ,  Ejector  OFF 
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Figure  H3  AO A  - 


,  Ejector  ON 


Figure  23  AGA 


Figure  26  AOA 


In  summary,  the  visual  effects  of  ejector  operation 

are : 

1 .  The  appearance  of  an  increase  in  pressure  differ¬ 
ential  due  to  the  upward  movement  of  the  streak  lines. 

2.  The  appearance  of  an  increase  in  effective  angle 
of  attack  and  increased  circulation  as  evidenced  by 
movement  of  the  leading  edge  stagnation  point. 

3.  The  appearance  of  a  more  well  defined,  less  turb¬ 
ulent  wake  region  above  the  aft  airfoil. 

1;.  The  evidence  of  increased  entrainment  due  to  the 
ejector's  employment. 

Based  on  an  analysis  of  the  flow  visualization  photo¬ 
graphs,  locations  for  obtaining  LDV  data  were  selected. 

The  LDV  results  for  AOA  =15°  and  AOA  =  0°  are  presented 
below. 

Laser  Velocimeter 
AOA  =  15°: 

The  first  LDV  results  presented  are  for  an  angle  of  at¬ 
tack  of  15°  with  the  ejector  duct  open.  Figure  28  shows  the 
survey  locations,  lettered  A  through  M,  which  were  conducted. 
Approximately  130  points  were  evaluated  both  with  and  without 
blowing.  Mean  velocity  profiles,  Um^  vs.  z,  and  turbulence 
intensity  profiles,  ETA  vs.  z,  are  presented  in  Figures  29 
through  54*  Please  note  carefully  the  z-ordinate  for  each 
profile:  z'  indicates  a  displacement  from  the  chord  line, 
while  z  indicates  a  displacement  from  the  surface  of  the  air¬ 
foil. 

The  first  survey  presented.  Survey  A  (Figure  29),  is 
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that  of  the  free  stream  conducted  at  a  point  0.65C  ahead 
of  the  airfoil's  leading  edge.  The  mean  velocities  were 
sufficiently  close  to  the  pitot  tube  established  tunnel 
speed  of  8  m/sec  to  give  confidence  in  the  pitot  tube  as 
a  means  of  setting  tunnel  speed.  Also,  to  within  the  ac¬ 
curacy  and  ability  of  the  LDV,  the  smoke  tunnel  free  stream 
exhibits  zero  turbulence  intensity  (Figure  42). 

In  reviewing  the  results  of  the  AOA  =  15°  Surveys  B 
through  M,  several  observations  can  be  made.  First,  the 
mean  velocities  with  ejector  operation  were  consistently 
greater  on  the  upper  surface  and  slower  on  the  lower  surface. 
Thus  a  greater  pressure  differential  exists  between  the  suc¬ 
tion  and  pressure  sides  of  the  overall  airfoil.  In  addition, 
the  higher  velocity  fluid  on  tne  upper  surface  is  better 
able  to  overcome  the  ad verse  pressure  gradient  and  thus  de¬ 
lay  separation.  Secondly,  with  ejector  employment  there 
was  an  increase  in  turbulence  near  the  airfoil , surf ace , 
coupled  with  a  more  rapid  decrease  in  turbulence  back  to 
the  free  stream  value  as  the  control  volume  was  displaced 
from  the  surface.  The  further  aft  of  the  leading  edge  the 
more  this  trend  intensified  as  evidenced  by  noting  tne 
turbulence  intensity  profiles  of  Surveys  F,  I,  and  K  (Fig¬ 
ures  47,  !?0,  and  52).  Also,  the  trend  is  evident  on  both 
the  upper  and  lower  surfaces,  as  Surveys  G  and  L  (Figures 
48  and  53)  will  illustrate.  This  increase  in  turbulence 
near  the  surface  with  ejector  operation  should  further 
energize  the  fluid  and  delay  separation.  The  observation 
of  a  more  rapid  decrease  in  turbulence  intensity  to  the 
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free  stream  value  with  ejector  operation  at  Surveys  K  and 
M  (Figures  52  and  54)  tends  to  confirm  the  flow  visualization 
interpretation  of  a  smaller,  more  well  defined  wake  region 
above  the  aft  airfoil.  Although  a  more  turbulent  fluid 
layer  suffers  more  viscous  drag,  a  delay  in  separation  coupled 
with  a  smaller  wake  region  should  improve  overall  performance. 
Also,  note  the  effect  of  the  converging  ejector  duct  on  the 
flow  near  the  surface  of  the  aft  airfoil  in  Surveys  K  and  M. 
The  higher  than  free  stream  velocities  indicates  greater  en¬ 
ergy  flow  on  the  aft  airfoil  surface.  The  effect  is  enhanced 
by  the  additional  momentum  flux  of  the  ejector.  A  third  ob¬ 
servation  at  this  AOA  concerns  Survey  J  (Figures  36  and  51 ) 
which  was  conducted  in  the  ejector  duct  downstream  of  tne 
ejector  face.  Without  the  ejector  on  an  increase  in  mean 
velocity  above  free  stream  is  noted  due  to  the  constricting 
duct  and  flow  over  the  aft  airfoil's  surface.  However,  with 
the  addition  of  the  ejector  flux  a  larger  mean  velocity  is 
obtained.  Additionally,  the  turbulence  level  is  greater 
with  ejector  operation.  This  indicates  an  increased  mix¬ 
ing  of  the  primary  ejector  fluid  and  the  secondary  entrained 
fluid.  Thus  higher  energy  flow  exits  the  duct  and  energizes 
the  aft  airfoil's  upper  surface  boundary  layer. 

Another  interesting  result  of  ejector  employment  is  il¬ 
lustrated  by  Survey  H  (Figure  }b) .  This  survey  was  conducted 
approximately  0.80  cm  aft  of  the  ejector  face  and  shows  the 
entrainment  of  secondary  fluid  into  the  ejector  duct  region 
by  the  ejector.  The  velocities  with  ejector  operation  are 
larger  in  the  areas  above  and  below  the  ejector.  AI30,  it 
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can  be  seen  that  without  the  ejector  on,  near  stagnation  con¬ 
ditions  are  present  at  the  ejector  face.  With  the  ejector  on, 
however,  the  established  ejector  velocity  of  approximately 
16.0  m/sec  is  attained.  Note  also  with  ejector  operation 
in  Figure  49  the  increase  in  turbulence  within  the  duct  and 
the  characteristic  more  rapid  decrease  in  turbulence  as  the 
control  volume  was  moved  away  from  the  airfoil. 

A  final  important  effect  of  ejector  operation  documented 
at  AOA  =15°  was  the  apparent  downward  or  counterclockwise 
movement  of  the  leading  edge  stagnation  point.  In  an  at¬ 
tempt  to  quantitatively  confirm  this  movement,  Surveys  D  and 
E  were  taken  (Figures  3 2  and  33) •  Surveys  D  and  E  are  loc¬ 
ated  1  cm  and  1 .4  cm,  respectively,  behind  the  leading  edge 
of  the  model.  Thus  Survey  E  is  located  further  down  the  air¬ 
foil's  surface,  i.e.,  in  a  counterclockwise  direction.  As 
previously  observed  in  the  flow  visualization  photographs  a 
trend  of  stagnation  point  movement  in  a  counterclockwise 
direction  is  noted.  Both  surveys  indicate  a  definite  de¬ 
crease  in  mean  velocity  with  ejector  operation.  Near  the  sur¬ 
face  at  Survey  D  the  velocity  without  blowing  is  near  zero, 
indicating  close  proximity  to  the  stagnation  point.  The 
mean  velocity  near  the  surface  with  ejector  operation,  how¬ 
ever,  is  negative.  Looking  at  Survey  E  the  velocity  near 
the  surface  with  blowing  is  still  negative  and  thus  has 
not  achieved  stagnation,  while  the  velocity  without  blow¬ 
ing  has  accelerated  past  the  stagnation  point  to  a  positive 
value.  This  can  best  be  visualized  as  in  Figure  27.  At  Sur¬ 
vey  D  the  velocity  without  blowing  is  nearly  zero,  whereas. 
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Figure  27  Leading  Edge  Stagnation  Point  Movement 

with  blowing  it  is  negative  and  thus  proceeding  against 
the  free  stream  apparently  on  its  way  to  flowing  over  the 
top  of  the  airfoil.  At  Survey  E,  however,  the  velocity 
without  blowing  has  accelerated  past  its  stagnation  point 
to  a  small  positive  value  on  its  way  to  flowing  beneath 
the  airfoil.  The  velocity  with  blowing  is  still  negative 
but  slightly  slower  as  it  accelerates  (although  in  a  neg¬ 
ative  direction)  to  Survey  D  and  then  o,rer  the  top  of  the 
airfoil.  Therefore,  it  can  be  observed  that  the  stagnation 
point  without  blowing  has  apparently  occurred  prior  to  Sur¬ 
vey  E,  while  the  stagnation  point  with  blowing  would  occur 
slightly  further  down  the  airfoil  surface  from  Survey  E. 
This  is  a  highly  simplified  model  of  the  actual  flow,  but 
the  trend  of  stagnation  point  movement  is  the  important 
point.  This  trend  indicates  a  stagnation  point  movement  in 
a  direction  that  would  yield  greater  circulation  with  ejec¬ 
tor  employment. 
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Figure  2b  Survey  Location*,  AOA 


igure  32  Mean  Velocity  Profile,  Survey  Location!  D  I  ,  AOA 


i^-ure  3k  Kean  Velocity  Profile,  Survey  Location  F  j  ,  AO  A 
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Figure  35  Mean  Velocity  Profile,  Survey  Location  I  G  ,  AOA 
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Figure  3?  Kean  Velocity  Profile,  ‘'urvev  ■.-cation 


40  Mean  Velocity  ?rof  Svrvsy  Location  I  I*  i,  AO  A 


Figure  1*1  Mean  Velocity  Profile,  Survey  Location 


Figure  42  Turbulence  Intensity  Profile,  Survey  Location 


Figure  43  Turbulence  Intensity  Profile,  Survey  Location 


Turbulence  Intensity  Profile,  Survey 


Turbulence  Intensity  Profile,  Survey  Location  H  I  ,  AOA 


Figure  50  Turbulence  Intensity  Profile,  Survey  Location 


Turbulence  Intensity  Profile,  Survey  Location  I  L  I  ,  AOA 


AOA  =  0°: 


The  next  angle  of  attack  to  be  investigated  with  the 
LDV  was  0°  .  Figure  56  shows  the  survey  locations,  let¬ 
tered  A  through  F  which  comprised  approximately  100  data 
points.  Since  the  relative  benefit  of  ejector  operation 
with  the  ejector  duct  open  was  documented  at  AOA  =  15°, 
the  investigation  at  AOA  =  0°  centered  on  comparing  the 
relative  performance  of  the  model  during  ejector  operation 
(ejector  duct  open)  with  an  equivalent  solid  airfoil. 
Therefore,  a  plexiglass  duct  blockage  element  was  designed 
and  fabricated.  With  this  blockage  element  and  black  tape 
the  model  could  be  temporarily  converted  to  a  solid  air¬ 
foil. 

In  an  attempt  to  determine  relative  drag,  wake  surveys 
(Survey  Location  F,  Figure  62)  where  conducted  as  follows: 

1.  Ejector  Duct  Open  (with  and  without  blowing): 

In  this  way  the  relative  drag  between  the  ejector 
on  and  ejector  off  operational  modes  with  the 
ejector  duct  open  can  be  compared. 

2.  Ejector  Duct  Blocked  (solid  airfoil,  no  blowing): 

With  this  survey  a  drag  comparison  can  be  made 
between  the  ejector  on  mode  with  the  ejector 
duct  open  and  an  equivalent  solid  airfoil. 

In  reviewing  the  data  presented  in  Survey  Locations 
A  through  E  (Figures  57  through  67),  which  compares  the 
model  with  ejector  duct  open  and  the  ejector  on  to  an  equi¬ 
valent  solid  airfoil,  the  same  performance  improvements  doc- 
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umented  at  AOA  =15°  can  be  observed: 

1.  Increased  pressure  differential  with  ejector 
operation  due  to  faster  mean  velocities  on  the 
upper  surface  and  slower  mean  velocities  on  the 
lower  surface. 

2.  An  apparently  smaller,  more  well  defined  wake 
region  above  the  aft  airfoil  with  ejector  operation. 

From  these  results  it  appears  that  the  ejector  wing  design 
is  superior  in  performance  to  an  equivalent  solid  airfoil. 

As  a  measure  of  the  relative  drag/thrust  augmentation 
of  the  ejector  wing,  wake  surveys  were  conducted  at  a  loc¬ 
ation  0.50C  downstream  of  the  trailing  edge  with  the  model 
at  0°  AOA.  The  surveys  are  presented  in  Figure  62. 

Wake  survey  data  can  be  used  to  find  the  drag  on  any 
body  through  a  control  volume  analysis  involving  the  equation 
of  motion  (Reference  15).  An  equation  for  the  drag  can  be 
written  as: 


Vmx)dZw  +  (mU) cos(9  ) 

v  © 


where  O' 


dzw 


(mU )  e 


Drag  per  unit  span  (or  thrust  per  unit 
span  if  negative) 

Air  density  =  0.0766  lb  /ft-* 

m 

Wake  survey  mean  velocity  component 

Coordinate  measured  normal  to  the  mean 
velocity  at  the  wake  survey  location 

Differential  area:  drw  x  (1),  i.e.,  unit 
span 

Momentum  flux  of  ejector 
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9  =  Angle  between  ejector  flux  and  the  free 

stream  with  the  model  at  0°  AOA,  0  =17° 

© 

Figure  55  illustrates  the  drag  equation. 


Figure  55  V.'ake  Survey 

The  first  term  in  the  drag  equation  represents  the  mom¬ 
entum  deficit  per  unit  span  between  station  1  and  station  2. 
The  second  term  represents  the  additional  momentum  flux  per 

unit  span  of  the  ejector.  Since  m  =  0.0U6  lb  /sec.  (mU)  = 

e  me 

O.315  lb^/ft.  Therefore,  with  the  ejector  off,  only  the  first 
term  is  used  to  calculate  the  drag.  But,  with  the  ejector  on, 
both  terms  are  employed.  In  thin  way,  the  additional  ejector 
momentum  is  accounted  for  In  the  control  volume  momentum 
balance.  A  drag  reduction  benefit  due  to  ejector  employment 
would  thus  be  evident  if: 

D'  ^  n« 

ejsctor  ON  ^  ejector  OFF 

Other  terms,  such  as,  pressure  differences  between 
station  1  and  station  2  due  to  the  tunnel  wall  friction  are 
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not  considered  significant  and  are  equally  present  in  both 
the  on  and  off  ejector  modes. 

Using  this  equation  and  graphical  integration  the  drag 
per  unit  span  of  the  model  can  be  calculated.  These  values, 
however,  are  relative  in  nature  and  are  used  here  only  for 
comparitive  purposes  to  show  the  effect  of  ejector  operation. 
Survey  Location  P  (Figure  62)  shows  the  effect  of  ejector  op¬ 
eration  with  the  ejector  duct  open  (ejector  ON  and  ejector 
OFF)  and  the  ejector  duct  blocked  (ejector  OFF).  Note  that 
the  velocity  deficit  is  decreased  with  ejector  employment. 

The  computed  relative  drag  values  are: 

D'  *  0.125  lbj./ft  (Ejector  OFF,  Duct  Blocked) 

D*  =  0.130  lbf/ft  (Ejector  OFF,  Duct  OPEN) 

D'  «  0.111*.  lbf/ft  (Ejector  ON,  Duct  OPEN) 

Thus  an  apparent  drag  reduction  for  ejector  employment  as 
compared  to  the  solid  airfoil  and  the  ejector  off,  duct  open 
operational  mode  is  evident. 
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Figure  56  Survey  Locations,  AO* 


igure  57  Mean  Velocity  Profile,  Survey  Location  A  I,  AO A 


Figure  61  Mean  Velocity  Profile,  Survey  Location  I  E  I,  AOA 


Turbulence  Intensity  Profile,  Survey 


Figure  66  Turbulence  Intensity  Profile,  Survey  Location 


A  Ejector  ON,  Duct  Open 


V.  CONCLUSIONS 


As  a  result  of  this  research,  conclusions  can  be 
drawn  as  to  the  aerodynamic  benefit  of  this  design. 

The  conclusions  that  relate  to  enhanced  lift  are  list¬ 
ed  first  followed  by  the  conclusions  that  relate  to  re¬ 
duced  drag.  Following  these  conclusions  are  three  add¬ 
itional  conclusions  that  are  related  to  both  increased 
lift  and  reduced  drag. 

The  conclusions  that  impact  lift  ares 

1 .  Ejector  employment  causes  an  increased  pres¬ 
sure  differential  between  the  upper  and  lower  sur¬ 
faces  of  the  airfoil. 

2.  The  effective  angle  of  attack  of  the  model  is 
increased  with  ejector  employment. 

3.  The  leading  edge  stagnation  point  appears  to  move 
down  the  forward  airfoil* s  lower  surface  with  ejector 
operation.  This  indicates  an  increase  in  circulation 
and  thus  an  increase  in  lift. 

The  conclusions  that  impact  drag  are: 

1 .  The  flow  on  the  upper  surface  of  both  the  forward 
and  aft  airfoils  is  increased  with  the  ejector  on. 
This  is  an  indication  that  the  ejector  energizes  the 
upper  surface  flows  and  thus  aids  in  fighting  adverse 
pressure  gradients.  This  may  delay  separation. 

2.  A  less  diffused,  more  will  defined,  and  smaller 
wake  region  appears  to  exist  above  the  aft  *.rfoll 
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as  a  result  of  ejector  employment. 


3.  With  ejector  operation  turbulence  levels  near  the 
surface  of  the  airfoil  are,  in  general,  increased. 

This  increase  in  turbulence  level  indicates  a  more 
energetic  boundary  layer  which  should  be  better  able 
to  overcome  adverse  pressure  gradients. 

4.  The  ejector's  momentum  flux  contributes  to  a  smal¬ 
ler  velocity  deficit  in  the  wake  region  (wake  survey 
location)  as  compared  to  an  equivalent  solid  airfoil. 

A  comparison  of  relative  drag  leads  to  the  conclusion 
that  the  ejector  provides  a  drag  reduction. 

The  conclusions  that  impact  lift  and  drag  are: 

1 .  The  ejector  entrains  secondary  fluid  flow  into  the 
ejector  duct  thus  increasing  the  mcraentum  flux  through 
the  duct.  This  should  further  enhance  the  airfoil's 
lift  increasing  and  drag  reducing  performance. 

2.  Flow  in  the  ejector  duct,  both  near  the  surface 
and  about  halfway  down  the  duct  is  more  turbulent 
with  ejector  operation.  This  increase  in  turbulence 
should  aid  mixing  of  secondary  and  primary  flows  and 
produce  a  more  energetic  flow  exiting  the  duct  and 
flowing  over  the  aft  airfoil. 

3.  The  ejector  wing  design  with  the  ejector  duct  open 
and  the  ejector  operating  appears  superior  in  performance 
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(increased  lift  and  decreased  drag)  to  an  equivalent 
solid  airfoil. 

On  the  basis  of  the  presented  results,  analysis,  and 
conclusions  it  appears  that  this  unique  ejector  wing  con¬ 
figuration  is  applicable  to  V/STOL  aircraft. 

Additional  conclusions  concerning  equipment  and  pro¬ 
cedures  are  as  follows: 

1.  Extreme  care  must  be  taken  when  aligning  and 
focusing  the  LDV  optics.  The  sensitivity  of  the 
optics  cannot  be  overemphasized. 

2.  Background  scattered  light  must  be  reduced  to 
the  maximum  extent  possible.  Steps  toward  tnis 
goal  include: 

a.  Insure  all  optical  surfaces  and  tunnel  walls 
are  clean  and  totally  grease  free. 

b.  Make  maximum  use  of  black  paper  and  cloth 
materials  to  reduce  reflections. 

c.  Insure  that  all  background  lights  are  ex¬ 
tinguished  prior  to  data  acquisition. 

d.  For  best  results  in  avoiding  scattered  laser 
light,  use  plate  glass  rather  than  plexiglass. 

3.  Due  to  the  limitation  of  1  MHz  maximum  phase  mod¬ 
ulation,  data  near  the  airfoil's  surface  (less  than  3 
mm)  and  in  other  regions  of  high  turbulence  (ETA>  60%) 
was  unobtainable. 
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VI.  RECOMMENDATIONS 


Based  on  this  study  which  utilized  a  LDV  and  the 
AFIT  smoke  tunnel  to  conduct  an  experimental  investi¬ 
gation  of  an  ejector  wing  design,  recommendations  for 
improvements  and  further  projects  are  presented. 

The  recommendations  are  divided  into  two  cat¬ 
egories:  those  relating  to  the  LDV  and  those  pertain¬ 

ing  to  the  ejector  wing  design  model. 

LDV 

1 .  For  ease  of  control  volume  translation  a  more 
sophisticated  setup  should  be  devised.  A  more 
stable,  larger,  and  most  importantly,  stationary 
primary  optical  bench  should  be  fabricated.  Upon 
this  stationary  bench  two  secondary  optical  benches 
could  be  used  to  achieve  translation  in  all  three 
dimensions.  The  incorporation  of  remotely  control¬ 
led  traversing  motors  would  greatly  aid  in  control 
volume  movement  and  should  help  alleviate  alignment 
and  maintenance  nightmares.  By  constructing  a  more 
stable  and  more  easily  aligned  optical  setup  the 
problem  of  floor/optical  bench  irregularities  which 
prevented  data  acquisition  very  near  the  airfoil's 
surface  may  be  avoided. 

2.  Incorporate  an  additional  laser  and  beamsplitter 
whose  control  volume  fringes  are  aligned  perpendicular 
to  the  original  fringes.  In  this  manner  two  components 
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of  velocity  can  be  determined  resulting  in  know¬ 
ledge  of  the  total  velocity  vector. 

3.  Increase  the  range  of  maximum  phase  modulation 
beyond  the  1  MHz  level.  In  this  way  data  in  areas 
of  high  turbulence  (ETA >60%)  could  be  obtained. 

4«  To  reduce  the  arduous  task  of  unskewing  the  dis¬ 
played  autocorrelation  function  and  subsequent  data 
reduction  process,  integration  of  the  LDV  system  with 
a  small  desk  top  computer  should  be  accomplished.  This 
project  is  being  conducted  in  parallel  with  this  in¬ 
vestigation  using  a  Hewlitt-Packard  (HP)  3052A  Auto¬ 
matic  Data  Acquisition  System.  This  project  should 
be  completed  as  soon  as  possible.  Once  the  entire  LDV 
system  is  automated  with  traversing  motors  controlled 
by  the  HP  for  control  volume  placement  and  data  reduc¬ 
tion  an  excellent,  real  time,  diagnostic  tool  will  be 
available  for  future  researchers. 

E.jector  Wing  Model 

1.  Incorporate  pressure  transducers  and/or  pressure 
taps  into  the  model.  In  this  way  pressure  distribution 
can  be  plotted. 

2.  The  model  is  equipped  with  trailing  edge  flaps  on 
each  airfoil.  The  forward  airfoil  has  a  12%  chord 

flap  and  the  aft  airfoil  has  a  19%  chord  flap.  The 
flaps  are  capable  of  +  4°  deflection.  Flow  field  var¬ 
iations  caused  by  the  deflection  of  these  flaps  could 


be  tie  basie  for  a  further  study  of  this  config¬ 
uration. 

3.  An  analytical  analysis  (finite  element,  panel¬ 
ing  or  other  computational  technique)  should  be 
conducted  on  this  configuration.  This  experimental 
study  could  be  used  as  a  guide  in  modeling  the  an¬ 
alytical  analysis.  In  this  way  the  experimental 
findings  and  the  analytical  results  could  be  com¬ 
pared  resulting  in  a  better  understanding  of  the 
ejector  wing  concept. 
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APPENDIX  A 


AUTOCORRELATION  FUNCTION  (Reference  4) 

An  autocorrelation  function  is  a  correlation  between  the 
same  fluctuating  random  data  measured  at  two  different  times 
at  the  same  point  in  the  flow.  Thus,  it  is  a  time  correlation 
rather  than  a  spacial  correlation.  As  an  example,  the  velocity 
at  a  point  in  the  flow  field  of  the  ejector  wing  may  vary  due 
to  turbulent  eddies  in  the  following  manner: 


Eac't  value  of  the  velocity,  u(t),  at  any  instant  in  time  is  the 
instantaneous  velocity.  An  estimate  of  the  autocorrelation  be¬ 
tween  values  of  u(t)  at  time  t  and  t  + X  is  obtained  by  taking 
the  product  of  u(t)  and  u(t  +  T)  and  averaging  over  an  observa¬ 
tion  time,  T.  Then  as  T  approaches  infinity  the  resulting  aver¬ 
age  product  will  approach  an  exact  autocorrelation  functir  de¬ 
fined  mathematically  as, 

T 

gu(T)  =  lim  £  \u( t )u( t+T)dt  =  u(t)u(t+T) 

T-*co  1  / 

0 
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In  order  to  nondimensionalize  the  autocorrelation  function, 
gu(T)  is  divided  by: 


which  ia  the  product  of  the  square  roots  of  the  time  averaged 
instantaneous  velocities  at  t  and  t  +  T . 


This  division  yields: 


G  (T) 


Since 


i,  i 

u(t)uU+T)  5s  [u2(t)]  [u2(t  +  X)] 


Schwartz's  Inequality 
(Reference  £1 ) 


then, 


-1  S  GU(T)  < 


A  typical  plot  of  the  autocorrelation  function  versus  delay 
time  (GuCT  )  vs  X )  is  shown  in  Figure  69  .  Note  that  Gu(0) 
has  been  normalized. 
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From  this  autocorrelation  function  mean  velocity  and 
turbulence  intensity  values  can  be  obtained.  However,  be¬ 
fore  proceeding  with  a  data  reduction  discussion  an  expla¬ 
nation  of  how  the  autocorrelation  function  displayed  on  the 
oscilloscope  is  transformed  into  a  useful  autocorrelation 
function  of  the  form  in  Figure  69  will  be  given. 

Since  the  oscilloscope  cannot  display  the  entire  auto¬ 
correlation  function,  ie,  x^oo  »  a  sample  of  the  function  is 
generated  by  operator  selection  of  an  appropriate  sample  time. 

In  general,  a  sample  time  is  selected  which  yields  the  best 
presentation  for  data  acquisition.  Sample  times  ranging  from 
0  to  9.99  seconds  in  .05  Jisec  increments  are  possible.  Select¬ 
ed  sample  times  are  composed  of  94  equal  time  steps  numbered 
5  to  99.  Thus,  tne  abscissa  of  the  oscilloscope  presentation 
represents  the  selected  sample  time.  This  sample  time  abscissa 
is  "free  floating"  on  the  display  requiring  the  operator  to  de¬ 
termine  its  location  and  orientation.  Location  and  orientation 
of  the  abscissa  will  be  discussed  in  detail  shortly.  The  or¬ 
dinate  is  the  number  of  photon  counts  detected  by  the  photo¬ 
multiplier  tube  and  processed  by  the  digital  photon  correlator. 
These  number  of  counts  (count  values)  on  the  ordinate  are  total¬ 
ly  arbitrary  and  will  vary  with  sample  time  and  oscilloscope/ 
correlator  settings.  The  arbitrary  nature  of  the  ordinate  is 
of  little  concern  since  the  function  will  be  eventually  nondim- 
ensionalized.  Also,  as  previously  stated,  the  sample  time  line 
(or  zero  ordinate  line)  is  not  displayed  on  the  scope  and;  there¬ 
fore,  must  be  determined.  This  determination  is  made  by  fitting 
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one  least  squares  second  degree  polynomial  curve  to  the  first 
three  maximum  count  values  and  a  second  least  squares  second 
degree  polynomial  curve  to  the  first  three  minimum  count  values 
as  shown  in  Figure  ?' 0  . 


Figure  70  Oscilloscope  Displayed  Autocorrelation 
Function 


The  arithmetic  mean  of  the  coefficients  is  then  calculated: 

A0+B0  Ai+E^  k2+b2 

Cq  =  — - —  ;  d  =  — - —  ;  C2  =  — - — 

A  new  mean  second  degree  polynomial  curve  is  then  drawn  as 
shown  and  this  represents  the  time  abscissa  (zero  ordinate) 
line.  That  is,  the  number  of  counts  corresponding  to  this 
line  is  considered  the  zero  count  line  and  all  other  count 
values  are  adjusted  accordingly.  This  time  abscissa  line  det¬ 
ermination  method  is  also  very  useful  for  correcting  a  skewed 
oscilloscope  display.  Skewedness  is  a  frequent  occurance  and 
is  due  to  either  an  insufficient  number  of  control  volume  fringes 
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or  detection  of  background  scattered  light.  Figure  71  shows 
a  skewed  oscilloscope  display. 


Number  of 
Counts 
( Arbitrary ) 


Pl(t)  =  A0+  Ait  +  A2t2 

Skewed  Autocorrelation  Function 
P3(t)  =  Cq+  Ci t  +  C2t2 


■P2(t)  =  B0+  Bi t  +  B2t2 


Sample  Time  (T) 

Figure  71  Displayed  Skewed  Autocorrelation  Function 


By  subracting  the  mean  curve  (P3)  from  the  minimum  (P2)  and 
maximum  (Pi)  curve  fits,  the  function  can  be  replotted  as 
anown  in  Figure  72  . 


Figure  72  Replotted  Unskewed  Autocorrelation 
Function 
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Once  an  unskewed  function  with  a  zero  ordinate  line  is  ob¬ 
tained  as  in  Figure  72  it  is  nondimensionalized.  This  is 
done  by  determining  the  difference  between  the  number  of 
counts  at  the  zero  ordinate  line  and  the  number  of  counts 
at  the  point  where  the  function  crosses  the  t=0  line  (time 
step  5)*  This  count  difference  is  then  divided  into  the 
count  values  at  all  points,  thus  nondimensionalizing  the 
entire  function.  The  resulting  function  is  the  desired 
autocorrelation  function,  Gu(t)  (Figure  73  )•  Note  that 
Gu(0)  has  been  normalized. 


Figure  73  Autocorrel a tion  Function 

A  computer  program  using  the  CDC  6600  computer  was  developed 
to  determine  the  Gu(t)  =  0  line,  correct  a  skewed  display  and 
assist  in  plotting  the  autocorrelation  function.  Once  the  os¬ 
cilloscope  displayed  autocorrelation  function  has  been  trans¬ 
formed  into  a  usable  function  data  reduction  is  the  next  step 
(Appendix  B  through  D). 


APPENDIX  B 


MEAN  VELOCITY  (Reference  6) 

By  definition  mean  velocity  is  the  average  velocity  at 
a  point,  with  respect  to  time. 

u  f  t ) 


t 

Figure  74  Mean  Velocity 


Mathematically,  the  mean  velocity  is  expressed  as: 


um 


=  u  =  1  i  m 

n—aj 


n 


i  =  1 


wlie re  n  is  the  number  of  time  samples  taken. 

In  order  to  determine  mean  velocity  from  the  autocorrelation 
function  some  preliminary  calculations  must  be  made.  First,  the 
fringe  spacing,  S,  must  be  determined. 


where  L  =  the  distance  from  the  laser  beam  intersection  point 

(control  volume)  to  a  perpendicular  reference  surface 

D  -  the  distance  the  two  beams  are  apart  at  the  reference 
surface 
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x  =  wavelength  of  helium-neon  laser  light  =  .6328  x  10"^  m 
J1  =  refractive  index  of  air  =  1.0 
Figure  75  illustrates  L,  D,  and  9. 


Then  the  formula  for  computing  the  mean  velocity  is: 


S 

UjnX  =  ( Peak-3  )(T) 


where  =  the  mean  velocity  component  parallel  to  the 

free  stream,  X 

Peak  =  the  time  step  for  the  first  peak  after  the 

first  valley  in  the  autocorrelation  function 
(Figure  7b  ) 

T  =  sample  time  selected 

Figure  7b  identifies  the  important  peaks  and  valleys  of  the 
autocorrelation  function  where: 

Gl  =  the  first  valley 

G2  =  the  first  peak  after  Gl 

G3  =  the  second  valley 
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Figure  76  Autocorrelation  Peaks  and  Valleys 


In  the  event  the  phase  modulator  had  to  be  employed  during 
data  acquisition  a  slightly  different  method  is  used  for  mean 
velocity  determination-  First,  the  included  half  angle,  9,  is 
calculated: 

bD 

tan  9  =  — 

L 

-1  JgD 
9  =  tan  — 

L 

Next,  the  mean  velocity  is  calculated  as  before: 


S 

(  Peak- 3 ) (T ) 


where  =  mean  velocity  with  phase  modulator  employment 

Then  calculate  the  doppler  shift  frequency  with  phase  modulator. 


2U,^(  sin9) 
_ 
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Next,  calculate  the  doppler  shift  frequency  without  the  phase 
modulator: 

Fd  =  F(i  1  A  F 

where  =  doppler  shift  frequency  without  phase  modulator 

employment 

A  F  =  the  frequency  shift  imposed  on  the  laser  beams 
by  the  phase  modulator 

The  sign  of  A  F depends  on  the  direction  the  fringes  are  made 
to  move  relative  to  the  flow.  If  they  move  with  the  flow.AF 
is  positive;  against  the  flowA  F  is  negative. 

Then  the  actual  mean  velocity  is: 

Fd  X 

UmX  2sinQ 

It  should  be  noted  at  this  point  for  completeness  tnat 
the  function  displayed  on  the  oscilloscope  can  be  used  to  det¬ 
ermine  mean  velocities  directly  witnout  correcting  for  skewed- 
ness  or  determining  a  zero  ordinate  line.  By  merely  determin¬ 
ing  the  time  step  at  which  the  first  peak  occurs  the  formula 
for  Um^  can  be  used  to  calculate  the  mean  velocity.  However, 
correcting  for  skewedness  and  determining  a  zero  ordinate  line 
is  necessary  for  the  calculation  of  turbulence  intensity  (Ap¬ 
pendix  C). 
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APPENDIX  C 


TURBULENCE  INTENSITY  (Reference  8) 

The  turbulence  intensity  parameter  is  a  means  of  meas¬ 
uring  the  velocity  fluctuations  from  the  mean  velocity  in  n 
turbulent  flow  (Figure  77  ). 


Figure  77  Velocity  Fluctuations 


Let  the  individual  fluctuations  about  um  be  represented  by: 


ui 

=  U1 

-  um 

=  u2 

um 

u3 

3 

li 

um  • 

e  tc . 

Then,  the  mean 

square 

of  u* 

*  Ums ,  is: 

_  n 

(u')2  =  lim  n  Y.  "  Ums 

n— "O  i=y 


The  root-mean-square  (RMS)  of  this  quantity  is  defined  as  the 


turbulence  intensity,  Urms: 


Turbulence  Intensity 


To  nondimensionalize  the  turbulence  intensity,  it  is  divided 
by  the  mean  velocity,  Um^  and  given  the  symbol,  ETA. 


ETA 


UR  MS 
Umx 


Tne  turbulence  intensity  is  obtained  from  the  autocorrelation 
function  by  the  following  method. 

First,  calculate  R: 


G2-G1  ,  .  n,  . 

R  =  G2TG3  (Figure  76  ) 

where  G1  =  the  first  valley 

G2  =  the  first  peak  after  G1 
G3  =  the  second  valley 

then. 


ETA 


%  [j1"- 


1  ) 


where 


■  -  ^  = 


number  of  fringes  in  one  half  of  the 
control  volume 


r q  =  laser  beam  radius  =  0.55  mm 


If  the  phase  modulator  is  employed  a  slight  variation  in 
the  above  equation  is  used: 
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ETA 


If  during  the  calculation  of  the  turbulence  intensity  an  un¬ 
defined  or  negative  value  is  obtained,  ETA  is  arbitrarily  set 
to  zero. 
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APPENDIX  D 


SAMPLE  CALCULATIONS 

In  order  to  illustrate  the  data  reduction  process  a 
sample  calculation  for  a  typical  survey  point  is  presented 
below. 

The  point  is  located  3*5  cm  above  the  airfoil  surface 
at  Survey  Location  K  with  the  model  at  15°  AOA.  Phase  mod¬ 
ulation  is  employed  and  a  skewed  oscilloscope  presentation 
was  obtained.  Tne  necessary  data  for  the  reduction  proce¬ 
dure  is: 

AOA  =  15° 

Ue/Uf3  -  2.0 

Sample  Time(T)  =  0.10  x  10”^  sec 

Phase  Modulation  =  201.9  KHz  (against  the  flow) 

Time  Step  Number  Digital  Count  Number  Peak(P)/Valley (V) 

5  32024546  p 

17  30973870  v 

31  31502864  p 

47  30982591  v 

56  31065118  p 

75  30885595  v 

Only  enough  data  to  allow  the  plotting  of  the  first  three 
autocorrelation  maxima  and  minima  are  needed  from  the  digital 
photon  correlator.  Therefore,  although  the  data  is  correlated 
and  stored  out  to  94  time  steps  (which  when  summed  together 
equal  the  sample  time  selected,  e.g .,  T  *  .10  x  10“^  sec)  only 


data  out  to  time  step  75  is  needed  in  this  case.  Also,  under 


certain  turbulent  conditions,  more  than  3  maxima  and  3 
minima  are  impossible  to  obtain.  Figure  ?8  is  presented 
as  an  illustration  of  this  point. 


Figure  78  High  Turbulence  Displayed 
Autocorrelation  Function 


The  first  step  in  the  data  reduction  process  is  to  un¬ 
skew  the  displayed  autocorrelation  function.  The  plot  of  the 
function  at  this  sample  point  is  shown  in  Figure  79  .  Next 
a  least  squares  second  degree  polynomial  curve  is  fit  through 
the  first  3  maximum  points  and  a  second  least  squares  second 
degree  polynomial  curve  is  fit  through  the  first  3  minimum 
points  as  illustrated  in  Figure  80  . 

Then  the  arithmetic  mean  of  the  coefficients  is  calcul¬ 
ated  : 

A0+B0 

Co  =  — ~ —  =  31  526669. 45 

Ai+Bi 

Ci  =  ~ —  =  -9123.61 
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Number  of  Counts  (x  10  )  (Arbitrary) 
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3.97 


A  2*^2 


Next  a  new  second  degree  polynomial  curve  is  then 
drawn  as  shown  in  Figure  8l  and  this  represents  the  zero 
ordinate  line.  The  skewed  displayed  function  can  now  be 
unskewed  by  subtracting  the  mean  curve  (P3)  from  the  minima 
(P2)  and  maxima  (P1 )  curve  fits,  and  replotting  as  illus¬ 
trated  in  Figure  62  . 

Finally,  to  complete  the  transformation  to  a  nondim- 
ensionalized  autocorrelation  function  the  entire  function 
is  divided  by  the  count  difference  at  Time  Step  5*  This 
value  is  543397-34  (Figure  83  ). 

With  the  displayed  function  unskewed,  mean  velocity  and 
turbulence  intensity  can  be  calculated  as  demonstrated  below 
(See  Appendices  A,B,  and  C  for  equation  details). 


Mean  Velocity 


S 


LX 

DU 


sin  0  =  tan  0 


L 


L  =  56.250  in 
D  =  1 .125  in 

X  =  0.6328  x  10’6  m 

U  =  1.0 


(56.25  in) ( .6326  x  10'6m) 
(1.125  in)  (1.0)  ” 


31.64  x  10“6  m 


sin  0 


(.5)0.125  in) 
(56.25  in) 


0.010 


105 


Number  of  Counts  (x  10  )  (Arbitra 


5  10  15  20  25  30  35  i|0  U5  50  55  00  65  70  75 


Time  Step 

Figure  8l  Least  Squares  Fit (Sample  Calculation) 


>6 


Count  Difference  (x  10~5) 


Figure  82  Unskewed  Autocorrelation  Function 
(Sample  Calculation) 
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J 


u. 


31 .61+  x  10  m 


mx 


Fd 


( Peak-3 ) (T ) 
2Ufjix(  sin  9) 


( 1 7-3) (0. 10  x  10‘6sec) 


=  11 


2(11.3  w/sec ) ( .01 )  _  ^£7  1 
.6328  x  10‘fem 


Pd  =  Fd 


U 


A  F 


mx 


=  357.1  KHz  -  201.9  KHz  =  155-2  KHz 

_  Fd  X  _  (155. 2/sec)( .6328  x  10~fem)  = 
2sin9  2 ( . 01 ) 


Turbulence  Intensity 

G2  -  G1  0.1+696  -  (-0.7339)  _  , 

G2  -  G3  0.4696  -  (-0.2283) 

»  .  £2  .  .$>  1  10~3".  =  ,7.383 
S  31  .61+  x  10“t>m 


— ~  =  0.00165 

2N2 


0.127 


.3  m/sec 


KHz 


m/sec 
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